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ABSTRACT 

A long-standing problem for models of galaxy formation has been the mismatch between 
the predicted shape of the mass function of dark matter halos and the observed shape of the 
luminosity function of galaxies. The number of massive halos is predicted to decrease as 
a power law (N oc M~ 2 ) out to very large masses, while the galaxy luminosity function 
cuts off exponentially at luminosities above L*. This implies that the efficiency with which 
gas cools onto massive systems is lower than expected. This letter investigates the role of 
radio-loud active galactic nuclei (AGN) in continually re-heating the cooling gas. By com- 
bining two observational results, the time-averaged energy output associated with recurrent 
radio source activity is determined, as a function of the black hole mass of the host galaxy: 
H = 10 214 (Mbh/-Wo) 1 ' 6 W. It is shown that for massive elliptical galaxies this radio-source 
heating balances the radiative energy losses from the hot gas surrounding the galaxy. The re- 
current radio-loud AGN activity may therefore provide a self-regulating feedback mechanism 
capable of controlling the rate of growth of galaxies. 
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1 INTRODUCTION 

In the widely accepted hierarchical clustering models of galaxy for- 
mation, the growth of structure is controlled by the gravitational 
collapse of haloes of dark matter, which merge together and accrete 
more mass to form progressively larger structures. These models 
have provided an excellent framework for describing the distribu- 
tion of mass in the Universe. The baryonic material, out of which 
stars and galaxies are formed, is superimposed upon this dark mat- 
ter distribution. As gas falls into the gravitational potential well of 
dark matter haloes, it is initially shock-heated to the virial temper- 
ature of the halo. This energy is then radiated away, and the gas 
cools and condenses, and begins to form stars. The properties of 
the resultant galaxies can be determined using semi-analytic mod- 
els (e.g. White & Frenk 1991); these adopt simple prescriptions for 
on-going physical processes, such as cooling of the gas, star for- 
mation, chemical enrichment, merging of galaxies, and feedback 
mechanisms, in order to predict the properties and distribution of 
galaxies. 

The galaxy luminosity function offers a powerful probe of the 
process of galaxy formation. In early, simple, semi-analytic mod- 
els, it was found that the predicted shape of the galaxy luminosity 
function differed greatly from that which is observed: many fewer 
galaxies are observed than were predicted at both high and low 
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masses. It was soon discovered that the problem of the lack of faint 
galaxies could be solved by including in the model the heating ef- 
fects of supernovae (e.g. White & Frenk 1991) and photoionisa- 
tion by the cosmic background (Efstathiou 1992), which both act 
to reduce star formation efficiency at the low mass end. The solu- 
tion to the deficit of bright galaxies was originally proposed to be 
due to cooling inefficiencies in massive galaxies, but these effects 
were subsequently shown to be insufficient. In order to reconcile 
the semi-analytic models with observations, some additional sup- 
pression of the gas cooling in massive galaxies was required. The 
generation of semi-analytic models of the late 1990's introduced 
this in an ad-hoc manner (e.g. by artificially switching off cooling 
in the most massive haloes), and were able to provide a very suc- 
cessful description of galaxy properties in the nearby and high red- 
shift Universe (e.g. Kauffmann et al. 1999, Somerville & Primack 
1999, Cole et al. 2000). In order to properly understand galaxy for- 
mation, however, it is clearly important to understand the physical 
processes involved in suppressing cooling in these galaxies. 

Active Galactic Nuclei (AGN) offer a promising solution to 
this problem (e.g. Benson et al. 2003, Croton et al. 2005, Bower 
et al. 2005, Scannapieco, Silk & Bouwens 2005). AGN are pre- 
dominantly found in the most massive galaxies (e.g. Kauffmann 
et al. 2003), and are sufficiently energetic that if their output could 
be efficiently coupled to the gas it would be more than sufficient to 
re-heat the gas and suppress further cooling. Prescriptions for AGN 
feedback have recently been introduced into semi-analytic models 
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of galaxy formation, and both Croton et al. (2005) and Bower et al. 
(2005) show that this feedback is able to solve the issue of the bright 
end of the luminosity function, whilst simultaneously solving other 
problems of galaxy formation models such as why the most mas- 
sive galaxies are so red. These prescriptions are necessarily simple 
in nature, however, and are not well grounded in observation. In- 
deed, the form of the AGN feedback adopted is very different in 
the two prescriptions (despite both providing good solutions) in- 
dicating that there is still much work to be done in understanding 
exactly how and when AGN feedback is important. 

The focus of this letter is to determine the AGN feedback due 
to radio-loud AGN. These are of particular interest because the ex- 
panding radio source provides a direct way for the AGN output to 
be coupled to its environment. Indeed, in clusters of galaxies con- 
taining powerful radio sources, X-ray observations have revealed 
bubbles and cavities in the hot intracluster medium, evacuated by 
the expanding radio source (e.g. Bohringer et al. 1993, Carilli et al. 
1994, McNamara et al. 2000, Fabian et al. 2003). The (pV) energy 
associated with these bubbles (which in turn arises from the me- 
chanical energy output of the radio jet) can be sufficient to balance 
the cooling losses of some clusters, at least for a short period of 
time (e.g. Fabian et al. 2003; Birzan et al. 2004). Radio-loud AGN 
have relatively short lifetimes (10 7 — 10 8 years), however, and so 
the balancing of cooling will not be long-lived. Nonetheless, Best 
et al. (2005a) recently showed that over 25% of the most massive 
galaxies show radio-loud AGN activity, which implies that AGN 
activity must be continually re-triggered. 

In this letter, recent observational results are combined to 
show that the time-averaged effect of heating by recurrent radio 
source activity balances the radiative energy losses in the hot en- 
velopes of the ellipticals, across a wide range of black hole masses. 
Section 2 provides a description of the observational results upon 
which this work is built. In Section 3 the results are discussed, 
the observations are compared with the prescriptions used in the 
semi-analytic models, and conclusions are drawn. Throughout this 
work, the cosmological parameters are assumed to have values of 
n m = 0.3, n A = 0.7, and H = 70kms- 1 Mpc~ 1 . 



2 METHODOLOGY 

2.1 The black-hole mass dependent radio luminosity function 

Best et al. (2005b) constructed a sample of 2215 radio-loud AGN 
with redshifts 0.03 < z < 0.3, from the second data release 
(DR2) of the Sloan Digital Sky Survey (SDSS), by comparing this 
with a combination of the National Radio Astronomy Observatory 
(NRAO) Very Large Array (VLA) Sky Survey (NVSS; Condon 
et al. 1998) and the Faint Images of the Radio Sky at Twenty cen- 
timetres (FIRST) survey (Becker et al. 1995). Best et al. (2005a) 
compared these AGN with the radio-quiet galaxies from the parent 
SDSS catalogue to derive the fraction of galaxies that are radio- 
loud as a function of black hole mass (Mbh, as determined from 
the host galaxy velocity dispersion; e.g. Tremaine et al. 2002) and 
radio luminosity separately. They found that the fraction of galax- 
ies that are radio-loud depends very strongly on black hole mass, 
going as / ra dio-ioud oc M^. On the other hand, the distribution 
of radio luminosities was found to be independent of black hole 
mass. Best et al. were able to parameterise this result using a mass- 
dependent broken power law model, such that overall the fraction 
of sources that are radio-loud AGN brighter than some 1.4 GHz 
radio luminosity L is given by: 
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Figure 1. The mechanical luminosity of radio sources, as estimated by 
Birzan et al. (2004) from cavities and bubbles in the X-ray gas, versus the 
monochromatic 1 .4 GHz radio luminosity. The solid line shows an error- 
weighted best-fit to the relation (in linear space). The errors on the me- 
chanical luminosities come from different assumptions in deriving the ages 
of the bubbles (see Birzan et al. for details). 
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with best-fit parameters of fo = (3.5 ± 0.4) x 10" 
0.1, (3 = 0.37 ± 0.03, 7 = 1.79 ± 0.14, L, = (3.2 ± 0.5) x 
10 24 WHz _1 . 

In a model in which radio source activity is continually re- 
triggered, this result can instead be interpreted (probabilistically) 
as the fraction of its time that a galaxy of given black hole mass 
spends as a radio source brighter than a given radio luminosity. 



2.2 Radio luminosity vs mechanical energy output 

Monochromatic radio luminosity does not provide a good indicator 
of the mechanical energy output of a radio source. Bicknell (1995) 
estimate that the kinetic energy output of a radio jet is typically a 
factor of 100-1000 higher than the total radio luminosity of a radio 
source. In order to use the above results to determine the time av- 
eraged heating output of radio-loud AGN, it is therefore necessary 
to derive a conversion between radio luminosity and mechanical 
energy output. 

Birzan et al. (2004) studied the cavities and bubbles that are 
produced in clusters and groups of galaxies due to interactions be- 
tween radio sources and the surrounding hot gas. They derived the 
(pV) energy associated with the cavities; combining this with an es- 
timate of the age of the cavities provides an estimate of the mechan- 
ical luminosity associated with the radio source. Figure 1 shows the 
mechanical luminosities derived by Birzan et al. (summed over all 
cavities associated with a given cluster or group) plotted against the 
monochromatic 1.4 GHz radio luminosity of the associated radio 
sources. As Birzan et al. point out, there is no single scaling factor 
between radio luminosity and mechanical luminosity. Nonetheless, 
it is possible to obtain a reasonable fit to the data using a simple 
power-law relation (shown as the solid line in Figure 1): 

t / T \ 0.40±0.13 

1 h = (3.0 ±0.2) ( _t^f-i ) ( 2 ) 
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where the errors are calculated using bootstrap analysis. 
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The scatter in Figure 1 is larger than the error bars, which 
suggests that there must be some additional dependencies in the 
conversion from radio to mechanical luminosity. This is not sur- 
prising since even for a source of fixed jet kinetic power the radio 
luminosity changes as the source ages (e.g. Kaiser et al. 1997). Im- 
portantly, however, the offset from the best-fit line in Figure 1 does 
not correlate with other properties of the host galaxy (optical lumi- 
nosity, velocity dispersion, environment) indicating that for all host 
galaxies, Eq. 2 provides a reliable measure of the mean mechanical 
luminosity associated with a source of given radio luminosity. 

Combining Eqs. 1 and 2 then gives the (probabilistic) frac- 
tion of time that a galaxy of given black hole mass spends pro- 
ducing a radio source of given mechanical luminosity, and from 
this the time-averaged mechanical energy output of a galaxy can 
be derived as a function of its black hole mass^ . This evaluates to 
H = 10 21 - 4 (M BH /M S ) 1 - 6 W. The black hole mass dependence 
of this heating rate arises entirely from the variation of the radio- 
loud fraction with black hole mass, whilst the radio luminosity to 
mechanical luminosity conversion determines the normalisation of 
the relation. 

2.3 Gas cooling rates in elliptical galaxy haloes 

The radiated energy losses from the haloes of hot gas surround- 
ing elliptical galaxies can be determined using X-ray observations. 
O' Sullivan et al. (2001) have compiled a (not complete, but rel- 
atively unbiased) catalogue of 401 early-type galaxies, for which 
they have tabulated the X-ray luminosities (Lx, in ergs^ 1 ), the 
optical B-band luminosities (Lb, in units of the solar luminosity 
in the B-band, Lbq = 5.2 x 10 32 erg s _1 ), and the T-type mor- 
phology parameter (de Vaucouleurs et al. 1976). The X-ray lu- 
minosities were derived using a common spectral model and dis- 
tance scale, and converted to a pseudo-bolometric band of 0.088 to 
17.25keV; O'Sullivan et al. show that extending this energy band 
further makes < 10% changes to the bolometric X-ray luminosity. 

The X-ray luminosity arises from two separate components 
of the galaxies: the sum of discrete sources within the galaxy (X- 
ray binaries, discrete stars, globular clusters, etc) and the cooling 
of the hot gas. The luminosity of the discrete sources, being stellar 
in origin, should scale roughly with the optical luminosity of the 
galaxy, and O'Sullivan et al. estimate that this scales roughly as 

« 10 29 ' 5 L B . This can dominate in low mass galaxies. For 
higher mass galaxies, the X-ray luminosity is dominated by the 
emission from the hot gaseous haloes, and O'Sullivan et al. find that 
Lx oc L\ (see also the review by Mathews & Brighenti 2003). 

For the current paper, analysis is restricted to those sources 

t Note that the conversion from radio luminosity to mechanical luminosity 
given in Eq. 2 is derived from radio sources in groups and clusters. It has 
been argued (e.g. Barthel & Arnaud 1996) that confinement of the radio 
lobes by the higher pressure intracluster medium leads to boosted radio lu- 
minosities compared to field radio sources of the same jet kinetic power. 
This might imply that Eq. 2 would underestimate the mechanical lumi- 
nosities for radio sources outwith clusters. However, the difference between 
cluster and field is only important for those sources whose radio emission 
extends beyond their host galaxy, and hence the extended surroundings con- 
trol the pressure of the radio lobes. In contrast, the radio sources which 
give rise to the bulk of radio source heating are low luminosity sources (cf. 
Fig. 3), which tend to be compact and more confined to the host galaxy. The 
radio lobe pressures are then more comparable to those of radio sources in 
clusters, meaning that Eq. 2 will be approximately valid. Note also that 
Eq. 2 gives values similar to those predicted by Bicknell (1995). 
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Figure 2. A plot of Lx versus L B for the elliptical galaxies in the sample 
of O'Sullivan et al. (2001). The large filled circles show the mean values 
of Lx for galaxies in 5 bins of Lg (each of width 0.4 in logLg); the 
mean values were calculated from both the observed values and upper lim- 
its, using the Kaplan-Meier estimator within the ASURV survival analysis 
package (LaValley et al. 1992). The dotted line gives the expected contribu- 
tion towards the X-ray luminosities from the sum of discrete X-ray sources 
within the galaxies. The solid line shows the predicted X-ray emission from 
the hot haloes of the ellipticals, assuming that heating from radio-loud AGN 
balances the cooling. 

with T < —4, that is, bona-fide elliptical galaxies. A plot of Lx 
versus L b for these galaxies is shown in Figure 2; the dotted line on 
that plot represents the average contribution to the X-ray luminosi- 
ties from the discrete sources. The residual effectively measures the 
cooling energy losses from the haloes of the ellipticals, as a func- 
tion of host galaxy luminosity. 

2.4 Black hole mass versus bulge luminosity relation 

In order to compare the average radio-AGN heating rates (as a 
function of black hole mass) with the average gas cooling rates (as 
a function of B-band luminosity), the black hole masses are con- 
verted to galaxy luminosities using the relation of McLure & Dun- 
lop (2002): log(M B H/M Q ) = -0.50M R - 2.96, where M R is the 
absolute R-band magnitude of the galaxy bulge. Because the sam- 
ple is restricted strictly to elliptical galaxies, the galaxy luminosity 
and the bulge luminosity are identical. The R-band magnitude is 
converted to a B-band luminosity using the average B-R rest-frame 
colour for elliptical galaxies, B-R«1.2 (e.g. Blanton et al. 2003), 
and the absolute magnitude of the sun (Mbq = 5.48). 



3 RESULTS AND DISCUSSION 

Figure 2 shows the Lx versus L b relation for elliptical galaxies 
from O'Sullivan et al. (2001). The solid line on this plot shows 
the time-averaged energy input into a galaxy by its recurrent radio 
source activity. The agreement between these heating and cooling 
rates across a wide range of host galaxy luminosities (masses) is 
remarkable. Averaged over time, the energy supplied to a galaxy 
by its recurrent radio source activity can balance the radiative en- 
ergy losses from its hot halo; thus, cooling of the gas is suppressed, 
and the radio source activity may control the rate of growth of the 
elliptical galaxy. 
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For this radio-source feedback model to work, it is neces- 
sary that the majority of the mechanical energy supplied by the 
radio activity is dissipated within the host galaxy halo. However, 
the most powerful radio sources can have linear sizes very much 
larger than their host galaxies. Combining the local radio luminos- 
ity function for AGN (e.g. Best et al. 2005b) with the conversion 
from radio to mechanical luminosity (Eq. 2), the contribution to 
the global mechanical heating rate from sources of different radio 
luminosities can be derived: this is shown in Figure 3. It is reas- 
suring that the heating is dominated by low luminosity sources, 
Li.4GHz 10 24 WHz _1 ; such low luminosity sources tend to be 
small, and their mechanical energy will indeed be confined predom- 
inantly to size-scales of the host galaxy halo. These low luminosity 
sources are also 'on' for a significant fraction of the time (up to 30% 
in the most massive galaxies; see Best et al. 2005a), meaning that 
the heat supply is reasonably continuous, as opposed to being in the 
occasional short-lived explosive events associated with the (rarer) 
more powerful sources. 

The volume-averaged mechanical energy heating rate due to 
all radio sources in the nearby Universe can be determined by in- 
tegrating Figure 3; because of the flatness of the relation at low 
radio luminosities, the total heating budget due to radio sources 
in the nearby Universe depends upon the radio luminosity down 
to which the radio luminosity function of AGN extends with- 
out change of sloped Adopting Li.4GHz = 10 22 WHz^ as the 
lower limit gives a zero-redshift volume-averaged heating rate of 
H = 4. Ox 10 31 WMpc~ 3 , but this value would nearly double if the 
slope of the radio luminosity function for AGN remains unchanged 
down to lO^WHz- 1 . 

Regardless of the precise value, the observed local heating rate 
due to radio-loud AGN is a factor of 10-20 lower than the local 
density of radio-mode heating in the Croton et al. (2005) semi- 
analytic model (9 x 10 32 WMpc~ 3 , as calculated from their Eq. 10 
and Fig. 3). Part of this difference is that the model cooling rates 
seem to be over-estimated by a factor of 2-3, due to the schematic 
nature of the underlying cooling flow model (White, Croton, private 
communication). The remainder of the difference arises because 
the level of radio heating in the Croton et al. model is tuned to bal- 
ance cooling in all massive galaxies, including the special case of 
brightest cluster galaxies at the centres of extreme cooling flows. 
The observations presented here, however, suggest only that radio 
source heating balances cooling radiation losses in 'typical' ellip- 
tical galaxies; there is no evidence for the additional heating that 
would be required to control cooling on cluster scales. An alterna- 
tive heating source is therefore required to solve the cluster cooling 
flow problem. This could still be associated with radio sources, but 
only if there is an additional mode of powerful radio source ac- 
tivity associated specifically with brightest cluster galaxies; this is 
beyond the test of the current observations. 

The balance between heating and cooling rates in typical el- 
liptical galaxies is likely to occur through a feedback-controlled 
episodic behaviour of the AGN activity (cf. Binney & Tabor 1995, 
Kaiser & Binney 2003): the radio source activity heats the sur- 
rounding gas; this then cools over time, until eventually some gas 



t Note that this result is very different from what would be obtained using 
the simple assumption that the mechanical energy scales linearly with the 
radio luminosity, when the bulk of the heating would arise from just a small 
range of source luminosities near the break of the radio luminosity function 
(cf. the dotted line in Fig. 3), and the range of radio luminosities considered 
would be less critical. 
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Figure 3. The global rate of heating in the nearby Universe produced by the 
sum total of all radio sources, as a function of their radio luminosity. The 
solid line shows the result obtained in this paper, deriving the mechanical 
luminosities of the sources from the fit to the data in Figure 1 . The dashed 
line shows the very different result that would be obtained assuming that the 
mechanical luminosity were directly proportional to the radio luminosity 
(the line is for L mcch = 300vL 1a q Hz , where the factor of 300 is roughly 
that estimated by Bicknell 1995). 

begins to condense down onto the central galaxy; a small propor- 
tion of that gas is driven down to the central regions of the galaxy, 
fuelling a new cycle of AGN activity; the resulting radio source 
re-heats the gas and cuts off the cooling once again. 

One interesting question is why, in this scenario, the radio- 
loud AGN fraction should scale as M^. Best et al. (2005a) found 
that the occurrence of low luminosity radio-loud AGN activity was 
completely independent of that of optical (emission-line) AGN ac- 
tivity, and therefore argued that radio-loud AGN were fuelled by 
the accretion of hot gas (as opposed to the cold gas which fuels 
optical AGN through a standard thin accretion disk). Croton et al. 
(2005) showed that the Bondi-Hoyle accretion rate of hot gas is 
"iBondi ~ G iim p kT 'Mbh / 'A (their Eq. 28), where \xm v is the 
mean mass of particles in the gas, and A is the cooling function. For 
an isothermal gas, T oc a 2 , which from the black hole mass versus 
velocity dispersion relation gives T oc M^. A is relatively inde- 
pendent of temperature (and hence black hole mass) at the temper- 
ature of elliptical galaxy haloes, and therefore the Bondi accretion 
rate scales roughly as Mg^, remarkably similar to the scaling of 
the radio-loud fraction: the higher proportion of radio-loud AGN 
in more massive galaxies may therefore arise naturally from the 
higher accretion rates. The physics underlying this accretion must 
be more complicated, since the Bondi accretion rate is for continual 
infall of gas whilst the mass-scaling of the radio-loud AGN frac- 
tion represents a scaling of the fraction of time that a given black 
hole is active, but it is certainly intriguing that the two scale with 
mass in the same way, and that this mass scaling matches that re- 
quired in the Croton et al. model to solve the current problems in 
the semi-analytic models. 

In conclusion, therefore, observations show that the time- 
averaged mechanical energy output associated with radio-source 
activity in massive galaxies is found to almost exactly balance 
the cooling radiation losses from the hot haloes of those galaxies, 
across a wide range of galaxy masses. This implies that the cooling 
of the gas may be feedback controlled by the radio activity. This 
feedback will control the growth rate of all galaxies in haloes in 
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which a quasi-static cooling halo of hot gas has been established; 
in practice this means all dark matter haloes more massive than 
2 - 3 x 10 n M Q (e.g. Bower et al. 2005, Croton et al. 2005). The 
mass-scaling of the mechanical energy output due to radio source 
activity matches both that predicted if they are fuelled by Bondi 
accretion of hot gas from the cooling haloes, and that required for 
radio-AGN feedback to successfully solve the problem of the over- 
growth of massive galaxies in semi-analytic models. 
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